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ABSTRACT 

The first luminous objects in the concordance cosmology form by molecular hydrogen cooling in 
dark matter dominated halos of masses ~ 10 6 Mq. We use Eulerian adaptive mesh refinement simula- 
tions to demonstrate that in the presence of a large soft ultraviolet radiation background, molecular 
hydrogen is the dominant coolant. Even for very large radiation backgrounds, the halo masses that 
cool and collapse are up to two orders of magnitude smaller than the halos that cool via atomic hy- 
drogen line cooling. The abundance of cooling halos and the cosmic mass fraction contained within 
them depends exponentially on this critical mass scale. Consequently, the majority of current mod- 
els of cosmological reionization, chemical evolution, supermassive black hole formation, and galaxy 
formation underestimate the number of star forming progenitors of a given system by orders of mag- 
nitude. At the highest redshifts, this disagreement is largest. We also show that even in the absence 
of residual electrons, collisional ionization in central shocks create a sufficient amount of electrons to 
form molecular hydrogen and cool the gas in halos of virial temperatures far below the atomic cooling 
limit. 

Subject headings: Cosmology: high-redshift — galaxy formation — star formation 



1. MOTIVATION 

Cosmic structure forms hierarchically. Any object 
in the universe today, started with copious numbers 
of small progenitors at redshifts currently inaccessi- 
ble to dir ect observations. Tr a ditionally in galax 
formation (Rces & Ostrikgi] [l977t IWhitc & Rccsl1l97i 



Dekel fc Reei Il987t IWhite fe Frenkl 119911: iBaugh et all 
20031) T„w = 10 4 K halos are assumed to be the 
first cooling halos. Nevertheless since the late 
1960's it has been known that molecular hydro- 
gen, formed in the gas phase, can dominate cool- 
ing in objects of smaller v irial temperatur e and 
mass (ISasiaw fe Zipovl H967L [Peebles fc Dickd U96l 
Yonevamal Il972t IHaiman et al.l 119961 : iTegmark et al.l 
1997t lAbel et alj 1 19981 I2000D Neglecting this early 



phase of H 2 cooling halos has been justified by arguing 
that H 2 is destroyed via radiative feedback effects (cf. 
Dekel fc Reedl 19871 IHaiman fc Loeblll997t IHaiman et all 
2000t iGlover fc Brandll2001l: iBromm fc Loebll2003D . The 
photo-dissociation of H 2 via the Solomon process by 
an early soft ultraviolet backgrou nd (UVB) is gener- 
ally assumed as the main reason (lOh fc Haimanl Bool 
ICiardi fc Ferrarall2005L IHaiman fc Brvanll2006f l ~ 

The mass scale of halos considered enters exponentially 
in the collapsed mass fraction and the abundance of ha- 
los. Figure Q] shows the predicted abundances of the ear- 
liest building blocks of galaxy formation as a function 
of redshift for the latest conc ordance cosmology using 
the Sheth-Tormen for malism (|Press fc Schechterl Il974t 
ISheth fc Tormenl 12002) . The different lines correspond 
to different virial masses. The solid line corresponds to 
halos with virial temperatures of 10 4 K, the temperature 
at and above which atomic hydrogen line cooling is domi- 
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nant. At redshift 30, e.g., the difference of abundances of 
2 x 10 5 M & and T vir = 10 4 K halos is five orders of mag- 
nitude. Even at redshift 10 this disparity is still a factor 
of a thousand. When studying reionization and chemical 
evolution of galaxies and the intergalactic medium, one 
needs to consider stellar feedback. The simple fact that 
the binding energy of the gas of smaller mass halos is even 
less than the kinetic energy deposited by even one super- 
nova (SN) is illustrated in Figure [TJ3. Surely whether 
the atomic hydrogen line (Lya) cooling halos are formed 
from pristine primordial gas or are mergers of many tens 
of progenitors that massive stars have enriched and ex- 
pelled the gas from should make a significant change in 
their further evolution. The minimum mass of star form- 
ing halos is undoubtedly an important issue independent 
of the techniques employed to study structure formation. 




Redshift 



Fig. 1. — Panel A: Sheth Tormen number density of dark 
matter halos as a function of redshift for T vir = 10 4 K, M = 
2 X 10 5 , 10 6 , and 6 X 10 6 M Q using WMAP 3 year data for param- 
eters. Panel B: Binding energies as a function of redshift for the 
corresponding halos (same line styles as in panel A) 
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Advances in cosm ologic a l hyd r odynamics a n d its 
numerical methods dCenl Il992t IZhang et all 119951: 
Katz fc Hernquisdll996tlAbel et al.lll997t lAnninos et all 
19971 iBrvan fc Normanl 119981 : iGnedin fc Abell 120011 : 
Ricotti et al.ll2002al )bl) allow now detailed investigations 
of all the relevant physical processes. Modeling the 
expected negative feedback from an early soft UVB is 
straightforward as a background flux only causes a spa- 
tially constant photo-dissociation rate in the chemical 
reaction network being solved when H2 does not exist at 
high e nough abundances to self-shield. Machacc k et al.l 
(|200lL MBA01 hereafter) used Eulerian adaptive mesh 
refinement (AMR) simulations to investigate the role of 
such a H2 dissociating (Lyman- Werner; LW) background 
on the minimum mass of halos within which primor- 
dial gas can first cool for a variety of radiation ampli- 
tudes. In addition to a LW background, the collapse 
of halos withi n relic H II reg i ons ca n be either delayed 
or catalyzed. M esigner et al.l (|2006l ) used AMR simula- 
tions with a short-lived 3 Myr hydrogen ionizing UVB 
that simulates a nearby massive, metal-free (Pop III) 
star. They found that halo collapses are prolonged if 
J21 2 ~ 0-1 and catalyzed if below this critical value, 
where Jff 2 is in units of 10 -21 erg s _1 cm -2 Hz -1 sr _1 at 
a wavelength of 912A. In the case of a large UVB, the 
collapse is delayed due to lower gas densities and higher 
cooling times. In the small UVB regime, excess free elec- 
trons in the relic H II region accelerate H2 formation. In 
both cases, feedback in relic H II subsides after ^30% of 
a Hubble time. Strong suppression of H2 formation also 
occur s in 10 6 MfT) halos wi th a LW background J^i > 
0.01. lYoshida et all (|2003l . YAHS03 hereafter) similarly 
addressed this issue using smoothed particle hydrody- 
namics (SPH). They found an additional effect on the 
minimum collapse mass of dynamical heating from the 
mass accretion history of the halo. As the heat input in- 
creases, the virial temperature must rise before H 2 cool- 
ing can start to dominate, and a cool phase develops in 
the center of the potential well. 

Self-consistent calculations in which the sources pro- 
duce the radiation backgrounds which in turn affect the 
number of new source s are feasible so far only with semi- 
analy tic approaches (jHaiman et alJ l2000t IWise fc Abell 
120051 WA05 hereafter) and small volu me cosmologi- 
cal sim ulations at low spatial resolutions (jRicotti et al.l 
2002a, b). From these studies, one can derive realistic 
upper limits on the amplitude of the expected soft UVB. 
In all studies that include radiation sources in halos less 
than 10 4 K halos, the largest the soft UVB flux can get 
before the T > 10 4 K halos dominate the emission is 
JW ~ 1 ( cf - iHaiman et all 120001: iRicotti et al1l2002alfbl . 
WA05). Interestingly, for a LW intensity of Jff ~ 0.1, 
MBA01 found that 2 x 1O 6 M halos were still able to cool 
and collapse. On the other hand at that J™ YAHS03 
suggest negative feedback should become so strong that 
the critical H2 fraction for cooling cannot be reached and 
cooling will not occur. However, they did not explore 
this further with detailed higher resolution simulations to 
check whether their analytical expectation would hold. 

We present a series of fourteen very high resolution Eu- 
lerian AMR simulations designed to see how the largest 
possible feedback may raise the minimum mass in which 
primordial gas will cool by molecular hydrogen. The sim- 



TABLE 1 
Simulation Properties 



Name H2 Residual e ^LW z a Zf, 



H2 Yes Yes 29.7 31.1 

H2LW22 .... Yes Yes 10" 22 28.3 27.5 

H2LW21 .... Yes Yes 1CT 21 24.4 24.7 

H2LW20 .... Yes Yes lO" 20 20.5 22.4 

noe-H2 Yes No 18.7 23.4 

noe-H2LW20 Yes No 1CT 20 16.8 21.4 

H+He No Yes 15.9 16.8 



Note. — These simulations are performed for both realizations. 



ulations techniques and details of the suite of calculations 
is the topic of the next section. In the following sections, 
we describe the results that show H2 cooling cannot be 
neglected in early structure formation. In the discussion, 
we describe the nature of the UVB and why H2 cooling 
can occur in such large radiation backgrounds. We also 
comment on the large range of questions in cosmological 
structure formation that this conclusion affects. 

2. SIMULATIONS AND ASSUMPTIONS 

We use the Eulerian AMR hydrodynamic code 
Enzo (jBrvan fc Normanl fl997. 1999) to study the impor- 
tance of H 2 cooling in early galaxy formation . Enzo uses 
an n -body adaptive particle-mesh solver (|Couchrnanl 
119911 ) to follow the dark matter (DM) dynamics. We 
perform two cosmological realizations with different box 
sizes and random phases and WMAP 1 year parameters 
of (h, A , n M , 0-8, n) = (0.72, 0.73, 0.27, 0.024ft- 2 , 
0.9, 1) (Sper gel et alj|2003f ). The significantly d ifferent 
third year WMAP (WMAP3: ISpergel et alj|2007l ) results 
favor lesser small-scale power that delays high-redshift 
structure formation by ~40% and alters th e statistical 
properties of DM halos (| Alvarez et al.ll2006bh . The ratio 
also only lowered by 5% to 5.70. However these 
differences have no effect on the evolution and assembly 
of individual halos studied here that have typical mass 
accretion histories. 

The initial conditions are the same as in IWise fc Abell 
(|2007f ). Both realizations have a top grid with a resolu- 
tion of 128 3 with three nested subgrids with twice finer 
resolution and are i nitialized at z = 12 9 (119 )^ with the 
COSMICS package (lBertschingeilll995l l200lh . The box 
size is 1.0 (1.5) comoving Mpc. The innermost grid has 
an effective resolution of 1024 3 with DM particle masses 
of 30 (101) M and a side length of 250 (300) comoving 
kpc. We refine the AMR grids when either the DM (gas) 
exceeds three times the mean DM (gas) density on the 
same level. We also refine so that the local Jeans length 
is resolved by at least 4 cells. 

We focus on the region containing the most massive 
halo in the simulation box and follow its evolution until 
it collapses to an overdensity of 10 7 that corresponds to 
a refinement level of 15 and a spatial resolution of ~3000 
(4000) proper AU. 

We perform each realization with seven sets of as- 
sumptions. Table [T] summarizes them. We use a 
nine-species (H, H + , He, He + , He + +, e~, H 2 , Hf , 
H~) non-equilibrium chemistry model (|Abel et al.lll997t 

t To simplify the discussion, simulation A will always be quoted 
first with the value from simulation B in parentheses. 
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Fig. 2. — Density-squared weighted projections in simulation A (left two columns) and B (right two columns) of the gas density (first 
and third columns) and temperature (second and fourth columns) at the times when the most massive halo starts to cool and collapse 
above an overdensity of 10 7 in the models. The rows show the H2, H2LW21, H2LW20, and noe-H2, noe-H2LW20, and H+He runs from 
top to bottom, respectively. Note the complex structure for the SimA-Noe-LW20 and SimB-Noe-H2 run in which central shocks lead to 
the formation of free electrons that promote the formation of H2 and triggering the collapse. The field of view in all panels is 1.2 proper 
kpc. The color maps are equal for all images. 
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lAnninos et al.l 1997) for all runs except the H+He runs 
that do not include H2 cooling. The nine-species runs 
are specified by "H2" a nd use the H 2 cooling rates 



are s pecified by t i^ a n 
from iGalli fe Pallal (fl998h . 



Above number densities of 
10 4 cm~ 3 or in an intense ultraviolet radiation field, the 
excited states of H2 become populated. The H2 colli- 
sional dissociation rates from lAbel et alj (|1997[ ) are cal- 
culated in the ground state; theref ore we use a density 
dependent H 2 dissociation rate from lMartin et all (1996) 
that considers this phenomenon. Runs with H2 disso- 
ciating (Lyman- Werner; LW) radiation are denoted by 
"LW" followed by its negative log-flux. We set Flw to 
1(T 22 , 1(T 21 , and 1CT 20 erg s _1 cm" 2 Hz" 1 because the 
first two are typical values one finds in semi-analytic 
models of reionization and the la tter investigates the 
case of an very large UVB (e.g. lHaiman et all l2000t 
iWise fc AbeJ 12005). We use the H 2 photo- dissociation 
rate c oefficient for the Solomon process from lAbel et al.1 
(fl997l) of fc diss = 1.1 x 10 8 F L w s . We do not consider 
the self-shielding of LW photons. Because the molecu- 
lar core only becomes optically thick in the late stages 
of collapse and above c olumn densities of 10 14 cm -2 
(jDraine fc Bertoldll 119961 ). we expect our results to not 
be drastically affected by neglecting LW self-shielding. 
Additionally, LW self-shielding may be unimportant up 
to column densities of 10 20 — 10 21 cm" 2 if the medium 
contains very large ve locity gradients and anisotropies 
(Gl over fc Brandll200l . 

Free electrons are necessary to form H 2 in the gas 
phase. In order to restrict H 2 formation to Lya line 
cooling halos in our "noe-" calculations, we reduce the 
residual free electro n fraction from ~ 10~ 4 (Peebles ! 19681 : 
IShapiro et al.l fl994"l to a physically low 10 -12 at the ini- 
tial redshift. This setup is designed to find the first halos 
that can collapse and form stars once free electrons from 
collisionally ionized hydrogen becomes ava ilable to cat- 
alyze H 2 formation (|Shapiro fc K ang 1983). 

This work is an extension of the original work of 
MBA01, adding the calculations with F L w = 1CT 20 and 
ones in which H 2 cannot cool until Lya cooling becomes 
efficient. We consider these extreme cases to strengthen 
the point made in MBA01 in which a UVB only increases 
the critical halo collapse mass, never completely sup- 
pressing the crucial importance of H 2 formation and cool- 
ing. Our maximum spatial resolution in the finest AMR 
level is a factor of four smaller than MBA01; however, 
this does not cause any differences between our work and 
MBA01 because these finest grid patches only exist in 
the dense, central core during the final 150 kyr of the 
collapse. 

2.1. Virial Temperature 

In galaxy formation models, the virial temperature is a 
key quantity as it controls the cooling and star formation 
rates in a given halo. We define a halo as the material 
contained in a sphere of radius r 2 oo enclosing an average 
DM overdensity A c of 200. For an isothermal singular 
sphere, the virial temperature 

_ ym p V L 2 
r ~ 2k ' 

where V} = GM/r-m o is the circular velocity (see 
iBrvan fc Normanll 19981 with 0=1). Here /i is the mean 
molecular weight in units of the proton mass m p , and k 
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Fig. 3. — Panel A: Halo masses of the most massive halos as 
function of redshift when they reach a central overdensity of 10 7 . 
This allows to translate the mass values in Panel B to be converted 
to cooling redshifts. It marks the runs with H+He (filled squares), 
H2 (open squares), -Flw = 10 -22 (open diamonds), -Flw = 10~ 21 
(filled diamonds), no residual electrons (open circles), and extreme 
feedback noe-LW20 (filled circles) runs. The two data points for 
each symbol represent simulations A and B. Even for the most 
extreme cases of feedback cooling occurs much earlier than in the 
atomic line cooling only case. Panel B: Central temperatures of 
most massive halo in the simulation as a function of its mass at 
different redshifts. The virial temperature computed form the dark 
matter halo mass at redshift 20 is the solid line. The dotted line 
is the fitted relationship between the central gas temperatures and 
the halo mass in models with residual electrons and H2 cooling. 



is Boltzmann's constant. We use this definition of T V i r in 
this paper with /i = 0.59. We choose this value of /x to be 
consistent with the literature on galaxy formation even 
though the halos presented in this paper are neutral and 
have n = 1.22. 

3. RESULTS 

We first describe the halo properties at collapse. Then 
we compare them to previous studies of collapsing halos 
in the presence of a soft UVB. 

3.1. Halo Properties 

Figure [5] shows density-squared weighted projections 
of gas density and temperature when each calculation 
can cool and collapse to an overdensity of 10 7 . It illus- 
trates the large difference in the sizes and morphologies of 
the collapsing halos in the various cases of negative feed- 
back. All panels have the same field of view of 1.2 proper 
kpc and same color scales. It is clear from the relative 
sizes of the collapsing halos that the critical halo mass 
to cool increases with the amount of negative feedback. 
The virial shock and numerous central shocks heat the 
gas to the virial temperature. The central shocks cre- 
ate fine structure seen in the temperature projections. 
In all of the H 2 cases, we see neither fragmentation nor 
large-scale disk formation. The internal structures of 
the halos with H 2 cooling and residual free electrons are 
similar to pr evious studies of Pop III star forming ha- 
los (MBA01: lAbel et alllmjol l200l iBromm et alJl200l 
lYoshida et al.l[2003l ). exhibiting a turbulent medium with 
a radially monotonically decreasing density profile and a 
cool central core. 
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Figure [3j^ depicts the halo mass and redshift when the 
halo collapses for all of the runs, and Figure [3j3 shows 
their central temperature at the same epoch. The col- 
lapse redshifts, z a and Zb, are also listed in Table Q] for 
simulations A and B, resepectively. As seen in other 
studies (MBA01, YAHS03), the minimum DM halo mass 
to collapse increases with the background intensity. The 
H+He case predictably collapses at T vil ~ 10 4 K, and 
all of the halos with H 2 cooling collapse at much smaller 
masses. The temperature of the central core increases 
with halo mass from 300 K to 1000 K for halo masses 
4 x 1O 5 M and 10 7 M Q . Restricting the data to models 
with residual electrons, the central temperature increases 
as a power-law, 

T c = AM vil B , (2) 

where A = 3.1±J;|, B = 0.355 ± 0.024, and M vir is in 
units of solar masses. This relationship is plotted in Fig- 
ure [3l 

With neither residual electrons nor an UVB (noe-H2), 
the most massive halo collapses at 9.8 (6.2) x 10 6 M Q 
at z = 18.7 (23.4). Here H 2 formation in the gas phase 
can only become important when sufficient free electrons 
are created by collisional ionization. Virial heating in 
the center of halos can increase temperatu res up to twice 
the virial temperature (Wise & Abel 2007) that collision- 
ally ionizes hydrogen in the centra l shocks and initiates 
H 2 cooling (|Shapiro fc Kandll987f ) in halos well below 
virial temperatures of 10 4 K. These shocks are abun- 
dant throughout the central regions. Figure 2] shows ra- 
dial profiles of temperature and electron fraction for both 
simulations and depicts gas shock-heating up to 2 x 10 4 K 
and raising electron fractions up to 10 -3 . The electron 
fractions remain at unrealistically low values less than 
10~ 6 in low density regions where gas has not been col- 
lisionally ionized. The higher density regions have con- 
densed to densities above 3x 10 2 cm -3 after free electrons 
in protogalactic shocks induced H 2 cooling. 

A similar but extreme model, noe-H2LW20, demon- 
strates that even in the presence of a very large UVB of 
Fun = 10~ 20 gas is able to form a cool and dense cen- 
tral molecular core at a mass of 2.7 (1.1) x 10 7 M Q at 
redshift 16.8 (21.4). Two major mergers in simulation 
A occur between z = 17-21, and the associated heating 
allows the halo to begin cooling by H 2 . A central core 
only forms once the system is adequately relaxed after 
the mergers, which causes the collapse mass difference 
between the realizations. 

By not fully resolving weak shocks in our main calcu- 
lations, it is possible to underestimate the electron frac- 
tion. We performed SimB-H2LW20 with an additional 
refinement criterion that resolves the "cooling length" , 
'cool = icooi/cs, by at least 2 cells. The large- and small- 
scale structure in the simulation is unchanged. When we 
resolve these weak shocks, the increased electron frac- 
tion marginally accelerates the collapse, which occurs 780 
kyr earlier at z = 22.5. The virial mass at this time is 
8.0 x 10 6 M Q compared with 8.4 x 10 6 M Q . Hence we be- 
lieve that the critical halo mass to collapse as a function 
of the LW background is independent of this refinement 
criterion. 

The combination of a recent major merger and colli- 
sional ionization produces complex structures as seen in 
the density and temperature projections of the SimA- 
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Fig. 4. — Radial profiles of temperature (top) and electron frac- 
tion (bottom) colored by density for the "noe-" simulations with no 
residual free electrons or UVB in simulation A (left) and B (right). 
The virial temperatures of these halos are 4600K and 4200K for 
simulation A and B, using equation JTJ with fi = 0.59. 

Noe-LW20 and SimB-Noe- calculations in Figure [U un- 
like the other H2 models with a single cool central core. 

3.2. Comparison to Previous Studies 

Through a series of AMR calculations with varying 
UVB intensities, MBA01 found the minimum DM halo 

mass 

M crit = 2.5 x 10 5 + 1.7 x IO^lw/IO" 21 ) - 47 M (3) 

in order to cool and condense 4% of the baryons. This 
fraction of cool and dense gas agrees with simu l ations 
of the formation o f Pop III stars (|Abel et al.l 120021 : 
lYoshida et aT1l2006D . There is some scatter of ~0.5 dex 
in this threshold mass (see also YAHS03). For the UVB 
intensities used in our models (-Flw = 0, 10~ 22 , 10 -21 , 
10~ 20 ), the critical collapse masses are 2.5 x 10 5 , 8.4 x 10 5 , 
2.0 x 10 6 , and 5.4 x 1O 6 M . Our calculations with 
H 2 cooling and residual free electrons agree with the re- 
sults of MBA01. 

YAHS03 studied the minimum collapse mass but also 
included the effects of self-shielding. Through their SPH 
simulations and arguments using equilibrium H 2 abun- 
dances, they conclude that an UVB intensity of J 2 i =0.1 
nearly prevents halo collapses below T V j r ~ 7000 K where 
Lya cooling becomes efficient. They also deduce that 
J 2i = 1.0 completely prevents any H 2 cooling in these 
low-mass halos, based on H 2 dissociation timescales. We 
find the contrary in our H2LW21 and H2LW20 calcu- 
lations where the most massive halo collapses with a 
mass of 4.5 (2.9) x 10 6 M Q and 8.4 (6.8) x 10 6 M Q , re- 
spectively. Even in our noe- runs, the halo collapses 
when T V j r ~ 4000 K, i.e. before Lya cooling becomes 
important, which is around the same mass scale that 
the H2LW20 runs condense. We ignore self-shielding in 
our calculations, but this would only decrease the critical 
collapse mass and strengthens our main conclusion that 
H 2 cooling is always dominant, even in the presence of a 
large LW flux. 

YAHS03 used cosmological SPH simulations and 
H 2 formation and dissociation timescales to argue that a 
LW background intensity of Jji > 0.1 suppresses H 2 for- 
mation so halos cannot cool before virial temperatures of 
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7000 K are reached. Employing the same argument, we 
see that the H2 formation timescale 



«H 2 



/h 2 



fc H -n H n e 0.92fc H -/ e n 



30 kyr, 



(4) 



with typical central values found in high-redshift ha- 
los before any rad iative cooling becomes efficient (see 
iWise fe Abellf2007h . Here / Ha = lO" 6 and /„ = 10~ 4 are 
the H 2 and electron number fraction, respectively, and 
n = 10 cm~ 3 is the baryon number density. /c H - « 10~ 15 
cm 3 s _1 is the H _ formation rate c oefficient by elec - 
tron photo-attachment at T = 1000 K (jAbel et al.lll997D . 
This timescale is a factor of 1000 smaller than the value 
calculated in YAHS03 because we use the quantities from 
the halo center as compared to the mean values. The 
H 2 dissociation timescale is kjA. = 23/J 2 i kyr, which is 
comparable with t^ 2 using the values above. 

The halo characteristics and the collapse redshift will 
likely depend on halo merger histories as seen in these 
two realizations. The better statistics of MB AO 1 sampled 
this effect well. Here the scatter of threshold mass is ~0.5 
dex and is smaller than the mass difference between halos 
with virial temperatures of 4000K and 10000K. Thus our 
limited sample of halos should not change our result of 
the importance of H 2 cooling in halos well below T v ; r = 
10 4 K, even with very large LW radiation backgrounds. 

4. DISCUSSION 

Structure formation in the high-redshift universe is 
contained within shallow potential wells that are sensi- 
tive to negative feedback from a UVB. Additionally local 
positive and negative feedback will influence star forma- 
tion and further complicate estimates of halo mass scales. 
Some examples include 

• Positive feedback — Enhanced H 2 formation in relic 
H II regions (e . g. lFerraral ll998t lO'Shea et aLll2005l; 
IJohnson et alJ l2007|) and ahead o f the H II ion- 
ization front (|Ricotti et alJ 120011: lAhn fe Shapirol 
20071) . dust and metal line coolingjfGloyerl 1200a 
Schneider et al.ll2006HJappsen et al.l[2007l ). 



• Negati ve feedback — Baryo n ic explusion from host 
halos (IWhalen et all 120041: iKitayama et "all 120041: 
lYoshida et all l2007t I Abel et all 120071) photo- 
evapor ation (|Susa fe U mcmura 2006) , entropy 
floors (IQh fe Haimanll2003D . 

These processes are not within the scope of this pa- 
per and will be considered in later publications that uti- 
lize three-dimensional radiation hydrodynamic simula- 
tions with Pop III star formation. Here we only focused 
on the effects of a UVB on low-mass halos. 

4.1. The Nature of the UVB 

The intensity of the UVB is a monotonically increasing 
function of redshift as more halos form stars. The UVB 
increases on the order of a Hubble time, which is much 
shorter than a dynamical time of a collapsing halo and 
justifies the use of a constant intensity in our calcula- 
tions. 

Self-consistent studies that evolve the UVB according 
to star formation rates only find J™ to be in the range of 
0.01 and 0.1 at redshifts 15-20 (YAHS03, WA05). WA05 



calibrated their model against the WMAP1 measurement 
of r = 0.17. With the WMAP3 result of the electron scat- 
tering optical depth r = 0.09 and less small-scale power, 
UVB intensities will be even lower at these redshifts. 

We can relate reionization to LW radiation by equating 
J21 in the LW band to a common quantity in reionization 
models, the ratio of emitted hydrogen ionizing photons 
to baryons, n 7j Hi/«-b, where ri& ~ 2 x 10~ 7 (1 + z) 3 cm~ 3 
is the cosmic mean of the baryon number density. As- 
suming that Jlw is constant in the LW band, the number 
density of LW photons is 



4?r 



Jlw 



dv 



= 1.19 x 10~ 5 J 2 L ^ cm" 3 , (5) 

where h p is Planck's constant and v\,v<x = 2.70 x 10 15 
Hz, 3.26 x 10 15 Hz bound the LW band. To relate J 21 to 
^7, hi/%, we must consider the intrinsic ionizing spec- 
trum and absorption from the IGM and host halo. At 
redshift 20, the majority of star forming halos host Pop 
III stars that emit a factor <^>hi — 10 more hydrogen ion- 
izing photons than LW photons because of its ~ 10 5 K 
surface temperature. Since the number density of sources 
exponentially increases with redshift, the majority of the 
early UVB at a given redshift originates from cosmolog- 
ically nearby (Az/z ~ 0.1) sources. Lyman line reso- 
nances absorb a fraction / a b s ~ 0.1 of the LW radiation 
in the intergalactic medium in this redshift ran ge, pro- 
ducing a sawtooth spectrum (Haiman et al. 1997). Addi- 
tionally, absorption in the host halo reduces the number 
of ionizing photons that escape into the IGM by a frac- 
tion f esc . For Pop H I halos, this facto r is close to unity 
(|Yoshida et alJ l2007t lAbel et all I2007D . By considering 
these multiplicative processes, we now estimate 



rc 7 ,m _ " 7 ,lw / 1 + z 
n b V 20 

fl + z 

= 0.64 J 2 i 



rib 



20 



0HI /esc / abs 



10 



(6) 



fahs 

"oT 



This estim ate is in agr e ement with the reionization 
models of lHaiman et alJ (|2000h and WA05 (see also 
iGnedin fe Ostrikeij 119971) . These models find that 
sources produce a large UVB of J 2 i ~ 1 prior to reion- 
ization. When Pop III stars dominate the UVB, the LW 
radiation will be small in comparison to the volume aver- 
aged hydrogen ionizing emissivity because of the intrin- 
sically hard Pop III spectra that peaks at ~300A. Hence 
high-redshift halos should not be exposed to a large UVB, 
i.e. J 2 i > 0.1, and H 2 formation will remain important 
before reionization. 

Nearby star formation can boost the LW radiation over 
its background value, but these bur sts are short-liv ed as 
Pop III lifetimes are only ~3 Myr (jSchaererl 12001 ). For 
example, a 100 M star produces 10 50 LW photons s _1 
and will produce > 0.1 in the surrounding 3 proper 
kpc, neglecting any H 2 self-shielding. 

The LW background is uniform outside these spheres 
of influence. The bursting nature of Pop III star forma- 
tion does not affect the time evolution of the background. 
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Fig. 5. — The most massive halo in SimB-LW20 (Flw 



10" 



10 10 

Radius [pc] 

5 ) twenty million years before the cooling core collapses. Slices of gas 



density (left) and the H2 formation timescale (Eq. [4]; right) through the densest point in the halo are in the top row. The bottom row 
contains a slice of temperature (left) and a radial profile of the H2 formation timescale, colored by gas density. The phase diagram and 
density slice have the same color scale. The slices have a field of view of 500 proper pc. In the central shocks, H2 formation timescales are 
lower than the dissociation timescale of 28 kyr with -Flw = 10 — 20 that is denoted by the horizontal line in the radial profile. The central 
core is now efficiently cooling and will collapse 20 Myr after these data. 



The intensity only depends on the number of sources in 
a redshift range Az/z = 13.6 eV / 11.18 eV - 1, where 
the two energies bound the LW band, because any radia- 
tion redward of the Lyman break contributes to the LW 
background. Using a conservative minimum halo mass 
for Pop III star forming halos of 3 x 10 6 M Q at redshift 
20, there are ~42000 halos that have hosted a Pop III 
star in the volume contained within Az, using WMAP3 
parameters with Shcth-Tormcn formalism. Clearly the 
background is uniform considering the sheer number of 
sources within this optically thin volume. Local pertur- 
bations from Pop III star formation should only affect 
the timing of nearby star formation but not the global 
star formation rate. 

4.2. H2 Cooling within a UVB 

Figure [5] shows SimB-LW20 twenty million years be- 
fore the central core condenses. At this time, the core 
is just beginning to cool by H2, catalyzed by the free 
electrons created in the central shocks. In these shocks, 
temperatures reach 1.4 x 10 4 K and electron fractions up 
to 1(U 3 exist there. These conditions result in H2 for- 
mation timescales less than 25 kyr, which is necessary to 
cool in a UVB of J21 ~ 1. Within the central 10 pc, hot 
and cold gas phases exist. The hot phase exists behind 
the shocks that have lower densities around 10 cm~ 3 and 
£h 2 < 25 kyr. This is where H2 cooling is catalyzed by 



collisional ionization in these shocks. The cold phase 
has already cooled through H2 and has high densities 
and larger <h 2 values. Both phases are apparent in the 
panels of Figure Similar conditions create H2 in the 
collapses in the "noe-" calculations, which have sufficient 
gravitational potential energy, resulting in temperatures 
above 10 4 K in central shocks. Hence H2 formation is 
possible in the centers of high-redshift halos with virial 
temperatures below 10 4 K, even with a UVB of inten- 
sity ~ 1, larger than expected from semi- analytic 
models of reionization. 

4.3. Impact on Semi-analytic Models 

Two consequences of a lower critical Lya cooling halo 
mass are more frequent and earlier galaxy formation and 
higher mass fractions in cooling halos. At redshift 20, 
e.g., abundances of T vir = 4000 K halos are an order 
of magnitude larger than T v - a = 10 4 K halos, result- 
ing from the exponential nature of Press-Schechter for- 
malism. The mass fraction contained in these halos is 
three times higher than 10 4 K halos. In semi-analytic 
models of reionization and chemical enrichment, the star 
formation rate (SFR) is linearly dependent on the col- 
lapsed mass fraction since the SFR is usually a product of 
mass fraction and star formation efficienc y, which is the 
fracti on of gas collapsing into stars fe.g. lHaiman et ail 
Il997f h The star formation efficiency for primordial stars 



WISE & ABEL 



is ~ 10 3 with a single massive star forming in dark 
matter halos with mass ~ 10 6 Mr, (lAbel et al.l 120021 : 
iBromm et all 120021: lYoshida et all |2006|). Thisfraction 
may rise to a few perce nt in dwarf galaxies as widespread 
star formation occurs (|Tavlor et alJlToM iGnedid [20001 : 
IWalter et alJ 1200 1). Various studies predict that a ma- 
jority o f the reionizing flux originates from dwarf galax- 
ies_(e.g. lCerj2003tlSokasian et alj2004HHaiman fc Bryan! 
200G). If the mass contained in star forming halos is three 
times greater than previously thought, some of the pre- 
dicted attributes, e.g. photon escape fractions and star 
formation efficiencies, of high-redshift dwarf galaxy will 
require appropriate adjustments to match observations, 
such as the WMA P3 measurement of optical depth to 
electron scatterin g dPage et al.l|2007f) and Gunn-Peter son 
troughs at z ~ 6 ([Becker et al.ll2001b iFan et al.H2002h . 

5. SUMMARY 

We conducted a suite of fourteen cosmology AMR sim- 
ulations that focus on the importance of H2 cooling with 
various degrees of negative feedback. We summarize the 
findings of each model below. 

1. The calculations with a UVB of F LW = (0, 10~ 22 , 
10~ 21 ) agree with the results of MBA01, where the crit- 
ical collapse halo mass increases as a function of UVB 
intensity. 

2. Above Flw = 10~ 21 , it had been argued that an 
H2 dissociating background would inhibit any H2 forma- 
tion until the halo could cool through Lya cooling. We 
showed that central shocks provide sufficient free elec- 
trons from collisional ionization to drive H2 formation 
faster than dissociation rates even in a Flw = 10~ 20 
background. 

3. In our "noe-" models, we explored when collisional 
ionization becomes important and conducive for H2 for- 
mation. This occurs at T vu ~ 4000 K. Recent major 
mergers above this mass scale create complex cooling 
structures, unlike the non-fragmented central cores in 
smaller halos. 

4. Even our most extreme assumptions of J21 = 1 
(.Flw — 10~ 20 ) and no residual free electrons cannot de- 
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feat the importance of H2 cooling in the early universe. 

O'S hea fc Normanl (|2007f ) independently studied halo 
collapses with Enzo and similarly considered primordial 
gas chemistry and nine different UVB intensities rang- 
ing from zero to J21 = 1. They agree with our conclu- 
sions in that primordial gas in T v ; r < lO^K halos can 
catastrophically cool and collapse even in models with 
J21 > 0.1. They attribute the collapse to the increased 
H2 cooling rates at higher temperatures that is caused by 
greater dynamical heating in halos with M v - iY > 10 7 Af Q . 
The cooling rate per molecule is 100 times larger at 2000 
K than at 500 K, typical of Pop III star-forming halos 
without an UVB. Most likely, the combination of the 
elevated H2 cooling rates and electron fractions from in- 
ternal protogalactic shocks instigate the halo collapses in 
a strong UVB (J 2 i > 0.1). 

In any case, H2 cooling triggers collapses in halos with 
virial temperatures well below 10 4 K. The lower critical 
halo mass, corresponding to T v ; r ~ 4000 K, increases 
mass fraction contained in these halos by three times at 
redshift 20 and the number density of high-redshift star 
forming halos by an order of magnitude! By considering 
additional cases of extremely large negative feedback, we 
have strengthened the results of MBA01 that H2 cooling 
plays a key role in high-redshift structure formation. We 
conclude that a UVB only delays and never completely 
suppresses H2 formation and cooling and subsequent star 
formation in these low-mass halos. 
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